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THE  object  of  this  paper  is  to  introduce  a  method  of  solving  simul- 
taneous numerical  functional  equations  which  is  original  and  is  thought 
to  be  new.  Any  process  of  solution  determines  the  value  of  each  variable 
as  the  sum  of  a  series  whose  terms  are  calculated  in  succession,  but  the 
practical  value  of  a  process  depends  not  merely  on  the  rate  at  which  the 
several  series  converge,  but  also  on  the  labour  involved  in  evaluating 
the  terms  ;  in  the  method  proposed  this  labour  is  almost  beyond  com- 
parison less  than  in  the  classical  method. 

The  problem  is  introduced  and  illustrated  by  a  particular  case,  arising 
out  of  a  showman's  puzzle  that  has  interested  a  large  number  of  mathe- 
maticians, and  when  the  general  method  has  been  explained  and  its 
application  to  two  special  problems  of  importance  indicated,  the  equations 
arising  from  the  puzzle  are  solved  and  some  interesting  consequences  of 
the  solution  described. 

I. 

A  familiar  figure  at  fairs  and  shows  is  the  sportsman  with  a  cloth  on 
which  a  large  circle  is  painted  and  five  smaller  equal  circular  discs  of  thin 
metal,  who  offers  the  holiday-maker  some  considerable  reward  if  he  can 
lay  the  five  discs  on  the  cloth  in  such  a  way  that  no  part  of  the  large 
circle  can  be  seen,  the  experimenter  of  course  paying  for  each  attempt. 
At  a  time  when  work  of  a  more  serious  kind  was  for  a  few  days  impossible 
to  me,  I  welcomed  the  problem  of  calculating  the  best  arrangement  of  the 
discs,  and  the  least  value  of  the  ratio  of  their  radius  to  that  of  the  painted 
circle,  and  the  results  may  interest  others.  The  equations  on  which  the 
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solution  of  the  problem  depends  were  much  easier  to  find,  and  much 
simpler  in  form,  than  I  had  expected,  but  the  numerical  resolution  of 
these  equations  with  sufficient  accuracy  seemed  likely  to  prove  intolerably 
tedious,  until  an  original  method  of  dealing  with  them  was  devised  : 
what  I  am  anxious  to  put  into  the  hands  of  mathematicians  is  this 
method. 

In  deciding  on  the  general  features  of  the  most  efficient  arrangement 
of  the  discs,  I  was  helped  by  having  a  specimen  of  the  apparatus  actually 
used  (belonging  to  Mr.  J.  H.  Grace,  whose  kindness  in  lending  it  both 
then,  and  on  the  occasion  of  the  reading  of  this  paper  to  the  Society,  I 
gladly  acknowledge)  ;  it  is  only  to  be  expected  that  no  great  margin  is 
left  to  the  inaccurate  speculator,  and  certain  types  of  arrangement  were 
seen  unmistakably  to  be  ineffective  ;  such  was,  for  example,  the  arrange- 
ment symmetrical  about  each  of  five  diameters,  the  small  circles  all  pass- 
ing through  the  centre  of  the  large  circle.  It  is  taken  for  granted  that 
there  is  symmetry  about  one  line,  a  common  diameter  of  the  large  circle 
and  of  one  of  the  small  circles.  If  K  is  the  centre  of  the  large  circle,  D 
the  centre  of  this  small  circle,  B  the  end  of  the  diameter  DK  of  the  large 
circle  which  is  not  covered  by  the  small  circle,  C  the  point  in  which  the 
small  circle  cuts  DB,  and  G,  H  the  points  in  which  the  small  circle  cuts 
the  large  circle,  the  arrangements  between  which  decision  must  be  made 
can  be  enumerated.  Two  circles  must  pass  through  B,  and  intersect  in  a 
point  L  in  DB,  which  may  be  identical  with  C,  or  may  be  a  distinct  point 
in  CB ;  let  one  of  these  circles  cut  the  arc  BG  of  the  large  circle  in  E,  the 
other  cut  the  arc  BH  in  F.  Of  the  remaining  circles  one  covers  E  and 
G,  the  other  covers  F  and  H.  If  L  is  distinct  from  C,  the  circle  covering 
E  and  G  covers  also  L  and  C,  and  either  passes  through  three  of  the  four 
points  E,  G,  L,  G  or  has  the  line  joining  two  of  them  for  a  diameter.  It 
L  coincides  with  C,  the  circle  BCE  cuts  the  circle  whose  centre  is  D  in  a 
point  M  distinct  from  C,  and  the  circle  covering  E  and  G  either  is  the 
circle  through  E,  G,  and  M,  or  has  one  of  the  lines  GM,  ME,  EG  for  its 
diameter.  It  would  be  possible  to  apply  calculation  to  each  case,  but 
actual  trial  is  sufficient  to  convince  that  the  only  arrangement  which  allows 
success  with  the  apparatus  used  is  of  the  last  type  ;  what  remains  for 
calculation  is  the  discovery  of  the  smallest  ratio  of  the  common  radius  of 
the  discs  to  the  radius  of  the  painted  circle  which  allows  this  most  effective 
arrangement  to  succeed,  and  the  determination  of  the  corresponding  posi- 
tion of  the  point  we  have  denoted  by  C. 

Let  b  be  the  radius  of  the  large  circle,  a  that  of  the  small  circles,  let 
c  be  the  distance  KD,  and  let  the  angles  BKE,  DKG,  BCE,  KDG  be  20, 
(f>,  \js,  TT  —  x-  It  is  a/b  that  it  is  our  chief  object  to  find,  but  (a  —  c)/6 


. »  * 
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derives  considerable  interest  from  being  extremely  small.     The  equations 
a  sin  \{s  =  b  sin  0,       (a—  c)sin  \fs  =  b  sin  (\fr—  20), 
c  =  b  cos  0  —  a  cos  x»       «•  sin  x  —  &  sin  (j>, 

implicitly  determine  6,  0,  \fs,  x  m  terms  of  c.  The  best  arrangement  of 
the  first  three  discs,  a,  b  being  supposed  constant,  is  that  in  which  c  has 
the  value  which  makes  EG,  FH  as  small  as  possible,  that  is,  is  found  by 
making  20  -|-0  a  maximum  subject  to  these  equations;  the  condition  for 


sin  2\/r  sin(x  —  0)  =  cos  x  !  1—  cos  20//-  —  0)J  , 

and  if,  when  the  position  is  such  that  the  five  equations  involving  c,  0,  <f>, 
^,  x  are  all  satisfied,  EG  and  FH  are  less  than  2«,  the  last  two  discs  can 
be  set  down  completely  to  cover  what  is  left  exposed  of  the  large  circle. 
The  least  value  of  a/b  which  allows  the  covering  is  that  in  which  EG,  FH, 
found  as  before,  are  equal  to  2a,  and  this  value  of  a/b  is  found  by  adding 
to  the  five  equations  already  written  down,  the  equation 

a  .  =  b  cos 


For  purposes  of  calculation,  the  ratios  a  :  b  :  c  are  eliminated  from  the 
equations,  and  the  resulting  equations  are  taken  in  a  form  involving  sums 
and  differences,  not  products  and  quotients,  of  circular  functions  ;  the 
equations  are 

^r-0  —  £0)  =  0, 


--  sin  (x~  0  —  £</>)  =0, 

2  sin  0  +  sin  (x+0)  —  sin  (x  —  0)  —  sin  ty  +  Q) 

—  sin  ty—  #)  —  2  sin  ty—  20)  =  0, 

cos  (2i/r  —  x  +  «/>)  —  cos(2i/r  +  x  —  <£)  —  2  cos  x 

+  cos(2^+x—  20)+cos(2x/.—  x—  20)  =  0, 

and  what  we  require  is  a  numerical  solution,  which  can  be  based  on  the 
fact  that  a  crude  approximation  is  given  by  the  values  36°,  36°,  72°,  72° 
for  0,  </>,  \/s,  x- 

It  is  evident  that  these  particular  equations  are  algebraic  and  even 
rational  in  the  tangents  of  the  angles  £0,  £0,  $\js,  $x<  so  that  an  algebraic 
equation  could  be  found  for  cos  (0  +  i0),  which  is  the  number  of  greatest 
interest  ;  to  calculate  by  Horner's,  or  any  other  of  the  familiar  methods, 
the  root  to  which  cos  54°  is  a  rough  approximation  is  in  theory  simplicity 
itself.  The  briefest  effort  will  convince  the  reader  that  it  is  impracticable 
to  solve  the  problem  on  these  lines.  The  alternative  is  to  regard  the 
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equations  as  functional,  rather  than  algebraic,  and  to  accept  suggestions 
from  the  differential  calculus  :  to  a  consideration  of  numerical  functional 
equations  we  now  turn. 

II. 

Of  all  the  common  problems  of  mathematics  whose  solutions  in  theory 
are  both  simple  and  complete,  probably  in  no  other  is  the  application  of 
the  theoretical  solution  to  a  numerical  case  as  tedious,  and  in  no  other 
does  this  application  find  the  mathematician  who  is  not  an  accountant  at 
as  serious  a  disadvantage,  as  in  the  problem  of  solving  a  set  of  numerical 
functional  equations.  Let  us  outline  the  classical  process  in  the  case 
of  three  equations 

F(x,  y,  z}  =  0,       G(x,  y,  z}  =  0,      H(x,  y,  z)  =  0, 
of  which  it  is  known  that  an  approximate  solution  is  given  by 

x  =  «ls       y  =  blt      z  =  Cj. 
If  the  corresponding  accurate  solution  is 

x  =  ctt+Xi,      y  =  &!  +  ?/!,      z  =  Cj 


an  approximation  to  the  values  of  xlt  ylt  zl  is  given  by 

*i  =  «4       IJi  =  b'2,       ZT_  =  Cz, 
where 

a'oF^a-L,  bv  c^-\-V^F,,(alt  blt  c^  +  CzF,(alt  blt  Cj)  =  —  F(alt  blt 


y(al,  blt  c^+dH^a^  blt  cj  =  —H(a1,  blt 
the  second  approximation  to  the  values  of  x,  y,  z  is  given  by 

x  =  r'2,      y  =  s'2,      z  =  4 
where 

+  a'2,       -so  =  b  +  b'z,       t'2  =  C  +  c^. 


To  obtain  a  third  approximation  to   the   values  of  x,  y,  z  the  process  is 
repeated,  terms  a's,  b'3,  c'3   being   obtained  from  three  equations  of   which 


the  first  is 


,     ,        = 


and  these  terms  being  added  to  those  already  known.     We  may  say  that 
the  exact  roots  r,  s,  t  are  the  sums  of  series  2«/TO,  S6,'t,  2c'n  whose  terms 
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are  calculated  in  successive  triads  from  three  equations  such  as 

f  i-CzVmt  Sm,  tnl)  —         r  (l',ni  Smi  t-init 


f'm,  s'm,  t'm  being  the  sums  to  m  terms  of  the  series  whose  w-th  terms  are 
a'm,  b'm,  c'm  and  the  first  terms  a[,  b{,  c{  being  the  given  first  approxima- 

tions «j,  &!,  Cj. 

This  process  is  open  to  a  twofold  criticism.  The  solution  of  simul- 
taneous linear  algebraic  equations  is  by  no  means  the  attractive  pro- 
cess in  practice  that  it  is  in  theory  :  if  there  are  only  a  few  equations 
and  the  coefficients  are  small  integers  the  solution  is  not  prohibi- 
tively troublesome,  especially  if  a  multiplication  table  or  a  mechanical 
multiplier  is  accessible,  but  given  coefficients  with  three  or  four  significant 
digits,  the  operation  is  excessively  tedious,  and  is  one  in  which  mistakes 
are  easy  to  make,  and  if  the  coefficients  are  irrational  it  is  troublesome  to 
decide  the  degree  of  accuracy  advisable  at  each  stage.  And  in  the  classical 
process,  although  the  coefficients  of  the  variables,  the  first  derivatives  of 
the  functions,  change  but  slightly  from  step  to  step  of  the  approximation, 
the  labour  of  solving  one  set  of  linear  equations  is  not  in  the  least 
diminished  in  virtue  of  the  work  done  in  solving  the  earlier  sets,  and,  of 
course,  increases  with  the  degree  of  accuracy  maintained.  In  the  alter- 
native process  which  is  to  be  described,  use  is  made  of  the  fact  that  the 
coefficients  vary  but  little,  and  subsequently  it  is  pointed  out  that  the  pro- 
cess is  applicable  if  the  coefficients  are  actually  constant  and  can  be  used 
in  the  solution  of  a  set  of  linear  algebraic  equations  with  effect  if  the  co- 
efficients are  complicated  or  irrational. 

If  the  coefficients  and  constant  terms  in  one  set  of  linear  algebraic- 
equations  differ  but  little  from  the  coefficients  and  constant  terms  in 
another  set,  the  values  of  the  variables  which  satisfy  the  one  set  differ 
from  the  values  which  satisfy  the  other  by  amounts  which  are  small  com- 
pared with  the  values  themselves,  it  being  assumed  that  the  determinants 
of  the  coefficients  are  not  small.  Thus,  if  \lt  /JL^  ...,  i>3  are  any  close 
approximations  to 

Fx(alt  h,  GI),    F,,(alt  blt  cj,     ...,    Hg(alt  bv  cj, 

the  values  of  the  nine  first  derivatives  of  three  functions  F(x,  y,  z),  G(x,  y,  z}, 
H(x,  y,  z)  for  values  av  b^  c^  of  the  variables  for  which  the  functions  are 
known  to  be  small  simultaneously,  and  if  a",  b",  c"2  are  such  as  to  satisfy 

the  equations  ,      »  .       7n  .       a  -,-,        , 

z  +  ViC*  --  —  F(al,  bit  erf, 


X2  a  "  +/x2  62  +  1/26-"  =  —  G(alf  blt  cj, 
S  =  —H(alt  blt  Cj), 
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then  aj+a'g,  &i~H&2i  ci~\~c'-2>  which  we  may  denote  by  r",  s",  t'2,  differ  from 
the  second  approximations  r'2,  s!2,  t'%  of  the  classical  solution  by  amounts 
of  the  order  g'[f",  where  /2"  is  the  greatest  of  the  moduli  of  the  three  terms 
a",  b'z,  c"  and  g"  is  an  approximation  factor*  depending  on  the  accuracy 
with  which  \,  /mlt  ...,  i/3  represent 

Fx(alt  blt  ca),    Fy(alt  blt  cj,     ...,    Hs(alt  blt  cj. 

Since  r'2,  s2,  t'2  differ  from  the  exact  solutions  r,  s,  t  by  amounts  of  order 
ft'2,  r'2,  s'2,  4'  differ  from  r,  s,  t  by  amounts  of  order  (</i'+/2")/^;  if/"  is  of 
a  higher  order  than  g",  the  approximation  furnished  b}7  r'!2,  s",  t"  is  not  as 
close  as  that  furnished  by  r'2,  s'2,  t'2,  but  the  somewhat  arbitrary  nature  of 
A1}  /uLlt  ...,  1/3  enables  us  to  take  for  them  rational  numbers  which  can  all 
be  comparatively  simple  without  the  value  of  g"  becoming  unduly  large, 
and  so  to  render  the  calculation  of  r",  s'2,  4  a  much  simpler  matter  than 
the  calculation  of  r'2,  82,  t'2. 

To  obtain  a  closer  approximation  than  that  given  by  r'2,  s'.2,  t'2,  we  have 
to  solve  linear  equations  which  either  are  equations  such  as 

tl  17>    I    'I         'I      J<\     I      7  "771     /     "          "       J."\     I        "  f    /     "          "       J."\  TV     "          "       4-"\ 

a3Fx(r2,  s.2,  t2)  +  b3Fv(rz,  sz,  t2)  +  c3Fz(r2,.s2,  tz)  =  —F(r.2,  s.2,  t2), 


or  are  equations  whose  coefficients  and  constant  terms  differ  little  from 
those  of  these  equations.  Now  it  is  assumed  throughout  that  none  of  the 
second  derivatives  of  the  functions  F,  G,  H  are  large  compared  with  the 
largest  of  the  first  derivatives,  and  it  follows  that  numbers  A1}  /m.^,  ...,  v3, 
which  are  approximations  to 

Fx(alt  blt  Ci),    Fy(alt  blt  Cj),     ...,    H,(al}  blt  cx), 
are  approximations  also  to 

Ft  •"    o"     t"\        7?  I  "      "    -f"\  TJ  I  •"      "    +"\ 

x\'2>  *2'    V2/'  2A    2>   ®2>    ^2/>       •••'       -favgi   *2'    "2'> 

and  indeed  to  the  values  of  the  first  derivatives  for  any  values  of  the  argu- 
ments not  differing  greatly  from  the  values  r,  s,  t,  which  we  are  endeavour- 
ing to  calculate.  Hence,  instead  of  solving  equations  such  as  that  last 
written,  we  may  solve  equations  such  as 

An 
ia3 


which   differ  only  in   the  constant  terms  from  the   equations   of  the  set 
solved  to  find  a'2,  b",  c.2;  the  sums  ?-3'»  s3>  4'>   ^na^  ig>  r2^~al> 


*  On  the  nature  of  this  factor,  see  Section  III  below. 
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differ  from  r,  ,s,  t  by  amounts  of  order  (g"  +/•")/$,  where  fl  is  the  greatest 
of  the  moduli  of  a3,  63,  c'3  and  g'l  is  an  approximation  factor  differing  but 
slightly  from  g". 

The  process  may  be  continued  indefinitely  :  the  solutions  are  found  as 
the  sums  of  series  2a"t,  26^,  Sc"ft  whose  terms  are  calculated  in  successive 
triads  from  sets  of  equations  such  as 

1,'t     II  'I          ill  \ 

^*"'  *m'     m'' 


the  arguments  r'n,  s'n,  t"m  being  the  ??i-th  partial  sums  of  the  series  2  a", 
26'^,  2c"i;  and  the  coefficients  Ax,  /x^  ...,  i/3  being  the  same  at  every  step  ; 
the  remainder  of  each  series  after  m  terms  is  of  order  (g'm-i~\~f%)fm>  where 
/,"  is  the  greatest  of  the  moduli  of  a"n,  b"m,  c"n,  and  g''n_l  depends  on  the 
accuracy  with  which  \,  /al,  ...,  i/3  represent 

777     /     II  II  ,11          \          7TT     /     It  II  ill  \  TT    /     II  II 

"xWm-U  S?/i-l'   tm-V'    -F  >J\rm-I>  Sm-l>   ^m-V'     •••>    -"aV^m-l'  Sm- 


m-V 


and  so  tends  as  m  increases  to  a  definite  limit  g"  dependent  on  the  near- 
ness of  the  coefficients  to  Fx(r,s,f),  Fy(r,s,t),  ...,  Hz(r,s,t).  Since/", 
tends  to  zero,  but  g"n_\  as  a  rule  does  not,  we  may  say  that  in  general  the 
remainder  of  each  series  is  of  order  </"/,"»  where  g"  is  a  fractional  approxi- 
mation factor,  though  in  exceptional  cases  the  remainder  may  be  of 
order  ff. 

When  we  have  to  solve  a  number  of  sets  of  linear  equations  with 
common  rational  coefficients  but  different  constants,  the  most  tedious  part 
of  the  work,  and,  what  is  equally  important  in  practice,  the  part  of  the 
work  in  which  the  greatest  care  is  needed  if  mistakes  are  to  be  avoided, 
can  be  performed  once  for  all.  From  the  set  of  coefficients  Al5  ^,  ...,  v.A 
a  set  />',',  a-'/,  ...,  r'3,  the  reciprocal  set  with  the  sign  of  every  member 
changed,  can  be  found  such  that  the  set  of  relations 


=  —  u,       ^a^       ^c  =  —  v,      ^ft--^       ^c  —  —  w, 

between  three  variables  a,  b,  c  and  three  variables  u,  v,  w  is  equivalent 
to  the  set  of  relations 

//          I         II         I         II  7  II          I         II          I         II  II          i         II         I         '/ 

a  =  plu-\-<Tlv-\-Tlw,    b  =  p.2u-\-(r2v-f-T2w,    c  =  p3u-\-<rsv-f-Tsiv, 
and  the  nine  coefficients  p'[,  <r">  ....T"    having  once  been  calculated,  the 
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successive  triads  of  which  we  are  in   search   are  obtained  from  equations 
such  as 


II  II        J."  \      I          "  /I  /     "  "         Jl  \     I          "  TT  /     "  "         J."  \ 

rm,  sm,  tJ+^Gfrn,  sm,  t'J+^H^,  sm,  tj. 

A  last  simplification  suggests  itself  immediately.  If  X1}  MI,  •••>  "a  are 
rational,  the  coefficients  p",  <r(,  ...,  T"  are  rational,  but  they  may  be  com- 
plicated fractions.  It  is  pointless  that  the  relations  such  as 

n      I      ii      i      '/ 
(I  =  pl  u  -|-  o-j  v  -f-  TI  w 

should  represent  the  relations  such  as 

\a-\-  jm1 


with  greater  accuracy  than  that  with  which  the  coefficients  \,  MI,  •  -.,  ^3 
represent  the  derivatives  Fx,  Fy,  ...,  Hz,  and  therefore  we  may  facilitate 
calculation  by  substituting  for  the  coefficients  p",  er",  ...,  -rl  any  coefficients 
plt  a-  1,  ...,  r3,  which  do  not  differ  greatly  from  them,  the  effect  being  to 
substitute  for  the  set  of  approximation  factors  g",  g'%,  ...,  a  set  of  approxi- 
mation factors  </!,  </2,  ...  tending  to  a  limit  g  which  may  be  either  larger 
or  smaller  than  the  limit  g".  So,  finally,  the  solutions  /•,  s,  t  are  the  sums 
of  series  2ani,  26,,M  2cm,  whose  terms  are  calculated  in  successive  triads 
from  the  equations 


)n,     m  --.  (T 


Sim   tm)  -f-  0"2  Gr  (/'TO,  S?n,, 

=  p3F(rm,  sm,  tm)-\-a-3G(rm,  sm,  tm)-\-r3H(rm,  sm,  tm), 


the  coefficients  pl}  a-^  ...,  r3  being  constant  throughout  and  simple  in 
form,  and  the  arguments  rm,  sm,  tm  being  the  m-th  partial  sums  of  the 
series  Sam,  2/jw,  2cni  themselves.  The  restriction  to  three  equations  in 
three  variables  has  been  purely  a  matter  of  convenience,  and,  in  general, 
we  have  the  theorem  : 

If  it  is   knoicn   that  an  approximate   simultaneous   solution   of  any 
number  n  of  independent  functional  equations 

Fp(xli  x2,  ...,  xn)  =  0     (p=  1,  2,  ...,  »), 
in  the  same  number  of  variables  is  given  by 

Xp  =  Opi     (p  =  1,  2,  ...,  M), 
rational  coefficients  pst  can  readily  be  obtained,  such  that,   if  sets  of 
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numbers  apm  are  calculated  in  succession  from  the  formula 

>i 

PqpFq(Tlm,,    )'-2,,i,     ...,    /',(,«)          (p    ~        1,    2,     ...,    /I), 


where 

»t 
/•jw  =    2  aik     (I  =  1,  '2,  ...,  n), 


fc=i 


00 


;t  series    2  a^  a/-e  convergent,  and  if  the  sum  of  the  series  2«y,,,t 

m=l 

is  ry),  Me  functional  equations  are  all  satisfied  for  the  set  of  values 
t\,  r2,  ...,  rn  of  the  variables;  the  series  are  ultimately  dominated  by 
geometric  series  with  a  common  ratio  dependent  on  the  differences  between 
the  values  chosen  for  the  coefficients  and  the  values  of  the  quotients 


( ?*        ?*n  V*    ^ 

unless  these  differences  all  happen  to  vanish,  in  which  case  the  dominating 
series  is  as  in  the  classical  method  of  approximation  a  series  of  the  form 


Theoretically,  the  classical  method  is  more  powerful  than  the  method 
described,  since  its  dominating  series  ultimately  converges  more  rapidly 
than  any  geometric  series,  but  in  practice  each  step  of  the  classical  method 
is  incomparably  more  troublesome  than  the  whole  group  of  steps  recom- 
mended here  for  advancing  the  approximation  by  the  same  amount. 

If  in  any  branch  of  applied  mathematics  the  solution  of  sets  of 
numerical  functional  equations  became  a  daily  necessity,  it  might  be 
worth  while  to  have  tabulated  sets  of  rational  coefficients  psf  correspond- 
ing approximately  to  standard  sets  of  coefficients  \t<  ;  the  undertaking 
would  be  weighty,  for  even  with  only  three  equations  the  number  of 
entries  would  be  very  large  if  the  number  of  different  values,  positive, 
zero,  and  negative,  which  each  coefficient  was  allowed  to  assume  was 
adequate. 

Before  returning  to  the  particular  set  of  functional  equations  for  whose 
solution  the  method  of  this  paper  was  devised,  let  us  refer  to  the  applica- 
tion of  the  method  to  the  two  simplest  cases ;  to  find  that  there  is  still 
something  to  be  said  on  familiar  problems  is  always  an  encouragement 
to  research. 
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III. 

We  have  already  remarked  that  if  the  coefficients  of  a  set  of  simul- 
taneous linear  equations  are  irrational  or  are  complicated,  the  solution,  to 
a  preassigned  degree  of  accuracy,  is  troublesome  :  we  can  apply  our 
general  method  to  deduce  the  solution  by  successive  steps  from  equations 
with  simple  coefficients.  Thus  to  solve  three  equations 


iz  =  P,      *x9*y*z  =  q,      a^^  =  r, 

we  take  any  numbers  X1}  /m^  ...,  »/3  which  do  not  differ  greatly  from 
/!,  <7ls  ...,  7?3,  and  tind  simple  rational  numbers  plt  v^  ...,  TS,  such  that 
the  equations 


between  two  sets  of  variables  x,  y,  z  and  u,  v,  10  are  approximately  equi- 
valent to  the  equations 

x  =  plu-\-a-lv-\-Tlw,     y  =  /o2w+<r2v-|-T2«0,     z  =  pBu-\-trsv-\-T3w  ; 
then,  if 


*  =  P—  (A 

q*  =  2  —  (/ 
r    =  r  - 


=  £2  — 


and  so  on,  the  series  #1  +  #2+...,  «/i  +  ?/2+  —  zv-\-z^-\-  ...  tend  to  the 
actual  values  of  x,  y,  z  satisfying  the  proposed  equations.  The  method  is 
particularly  useful  if  the  original  equations  have  irrational  coefficients, 
since  it  renders  it  unnecessary  to  determine  in  advance  how  far  accuracy 
in  approximating  to  the  values  of  the  coefficients  is  significant  in  the 
attainment  of  any  required  degree  of  accuracy  in  the  solutions,  and  it 
renders  it  possible  to  make  use  of  any  given  approximate  solution  in 
searching  for  an  approximation  still  closer.  So  true  is  it  that  when  co- 
efficients are  complicated  an  economy  is  effected  by  this  process,  that  the 
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simplest  way  of  conducting  an  approximation  to  a  solution  of  a  set  of 
functional  equations  by  the  classical  method  is  to  solve  hy  this  method 
the  various  sets  of  linear  equations  that  arise  :  this  mixture  of  methods  is 
naturally  less  satisfactory  than  a  frank  desertion  of  the  classical  method, 
but  when  in  any  case  a  point  is  reached  from  which  linear  equations  would 
complete  the  solution  to  the  required  degree  of  accuracy,  to  write  down 
these  linear  equations  enables  the  calculator  to  avoid  further  reference  to 
tables. 

In  the  case  of  linear  algebraic  equations  it  is  easy  to  shew  the  validity 
of  the  process  employed.     With  the  three  equations  just  discussed,  we  have 


\  1  +  (Pl/l  +  O'l/j  +  Ti/j)  ;   Zm+fcl^l  +  O-l^  +  T^g)  ,J  ,A+  (Pi  7^  +  ^  //o+M  7l3)  Z 

and  similar  expressions  for  ijm+i,  zin  +  \.     The  nine  coefficients 


can  be  made  as  small  as  necessary  by  proper  choice  of  the  nine  numbers 
Pn  o-p  ..,  T;J,  and  it  is  evident  that  if  the  greatest  of  their  moduli  is  not 
greater  than  \g,  and  g  is  a  proper  fraction,  the  series 


all  converge  not  less  rapidly  than  a  geometric  series  with  ratio  g.  In 
practice  the  method  may  often  be  found  to  succeed  even  if  the  greatest 
modulus  is  not  less  than  ^,  for  the  signs  of  the  various  terms  do  not  as  a 
rule  combine  in  the  most  unfavourable  manner  conceivable.  The  method 
being  proved  valid  for  linear  algebraic  equations,  its  validity  for  equations 
of  any  form  is  a  consequence  of  the  validity  of  the  classical  process. 

The  other  simple  case  of  which  we  wish  to  speak  is  that  of  a  single 
functional  equation.  What  we  say  is  that  if  %  is  a  first  approximation  to 
a  root  of  a  functional  equation 

F(x)  =  0, 

and  —p  is  any  number  not  very  different  from  the  value  of  the  reciprocal 
of  dF/dx  when  x  is  equal  to  als  then  the  root  itself  is  the  sum  of  the  series 
whose  partial  sums  are  calculated  in  succession  from  the  formulae 

«2  =  pF(al) }      as  =  pF(>:2) )      «4  =  pF(r:J        a5  =  pF(r4)  \ 


and  so  on. 
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For  example,  taking  2— #2  for  F(x)  and  1£  as  a  first  approximation, 
we  naturally  take  p  to  be  ^,  and  we  have 

a2  =  —    '083,     03  =  —    '0026,     «4  =  —    "00018,     %  =  —    '000006, 
r2  =       1-417,      ra  =       1-4144,      r4  =       1'41422,     r5  =       1'414214. 

*The  process  of  solving  a  single  functional  equation  may  be  illustrated 
graphically,  the  relation  of  the  method  of  this  paper  to  the  classical 
method  being  made  apparent.  The  problem  is,  given  a  point  A^  on  the 
axis  of  x  known  to  be  not  far  from  the  point  P  in  which  a  curve  y  =  F(x) 


meets  that  axis,  to  construct  a  sequence  of  points  Nlt  N%,  Na,  ...,  having 
P  for  limiting  point.  In  the  older  method,  from  N±  is  drawn  the  ordinate 
NlPl  to  the  curve  (the  line  parallel  to  the  axis  of  ?/,  not  necessarily  the 
line  perpendicular  to  the  axis  of  x),  and  the  tangent  Pi-/V"2  to  the  curve  at 
Pl  cuts  the  axis  of  x  in  the  second  point  N%  of  the  sequence  ;  A^  lies  on 
the  tangent  at  the  point  P2  in  which  the  line  Ar2P2  parallel  to  A^Pj  cuts 
the  curve,  and  so  on.  In  the  present  method,  Pv  is  found  from  Nlt  P.2 
horn  N.2,  and  so  on,  as  before,  but  the  lines  Pi-/V2,  Pa-^a,  •••  by  which  N>2 
is  found  from  P1?  N3  from  P2,  and  so  on,  are  not  tangents  to  the  curve,  but 
are  lines  all  parallel  to  some  direction  not  differing  greatly  from  that  of 
the  tangent  at  Pv  If  actual  tangents  are  used,  tables  must  be  consulted 
for  the  determination  not  only  of  the  lengths  of  the  ordinates  A^P^  -/V2P2,  . . ., 
but  also  of  the  slopes  of  the  tangents  P1N^  P.2Ar3,  ...  ;  if  the  lines  PiA^, 
P.2Ar3,  ...  have  a  common  direction,  the  triangles  Ar2Ar1P1,  N3N,2P2,  ...  are 
all  similar,  and  when  the  ratio  of  N^Ni  to  N^Pi  has  been  found  it  is  only 
for  the  values  of  the  ordinates  NiP^  N^PZ,  ...  that  reference  need  be 
made  to  tables.  This  graphical  consideration  brings  out  clearly  another 
point.  The  rapidity  of  the  approximation  depends  on  the  closeness  of  the 


*  The  remainder  of  this  section  has  been  added  since  the  reading  of  the  paper ;  the 
graphical  illustration  was  suggested  by  one  of  the  referees. 
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common  direction  of  PiN2,  P.2N.A,  ...  to  the  direction  of  the  tangent  at  P, 
the  point  to  be  found,  but  the  approximation  may  be  valid  even  if  the  two 
directions  differ  considerably  :  if  PPl  is  actually  a  straight  line,  and  N.2  is 


any  point  between  Nl  and  the  image  of  N1  in  P,  then  by  drawing  N>2  P2 
parallel  to  N^,  P,2N3  parallel  to  P^N.2i  N3P3  parallel  to  N^,  PSN+ 
parallel  to  Pi-/V2,  and  so  on,  we  can  construct  a  sequence  of  points  con- 
verging to  P,  the  convergence  being  the  faster  the  nearer  JV2  is  taken  to 
P.  Similarly  in  the  more  general  case,  there  may  be  much  latitude  in 
the  choice  of  the  common  direction  of  P\NZ,  P^NA,  and  the  rest,  and 
because  of  this  latitude  a  choice  may  be  made  which  gives  a  simple  value 
to  the  ratio  of  N^Nl  to  ^i-Pi»  °f  N3N.2  to  N.2P2,  and  so  on. 

Thus  we  see  a  connection  between  the  method  of  this  paper  and  a 
method  used*  in  a  number  of  problems  in  applied  mathematics.  The 
simplest  ratio  which  N.2Ni  can  bear  to  N1P1  is  unity,  and  to  draw  PiN.2  in 
such  a  direction  as  to  make  N.2Nl  equal  to  N1P1  is  to  evaluate  the  root  of 
the  equation  F(x)  =  0  by  the  steps 


a.2  =  — 
that  is, 

/•2  =  /-!  —     d-j),     r.A  —  r.2  — 

This  process  is  specially  useful   when  the  equation  to  be  solved  has  the 

particular  form  ,.,  , 

x  —  j  (x), 

for  if  we  write 

F(x)  =  x-/ 


we  see  that  the  approximation  then  advances  by  the  steps 
''2  =  /  O'i)  »     '3  =  f(>'J  »     ''i  —  /(*"£  > 


*  I  am  indebted  to  Dr.  Bi'omwich  for  drawing  my  attention  to  this  method, 
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which  are  so  simple  as  to  justify  considerable  slowness  in  the  convergence. 
It  can  be  shewn  that  convergence  is  certain  (though  it  may  be  slow),  if 
f'(x)  is  positive  and  less  than  unity  in  the  neighbourhood  of  the  root,  and 
since  if/"1^)  is  the  function  inverse  to  f(x),  the  equation  x  =  f(x)  is 
identical  with  the  equation  x=/~l(x)  and  df~l(x)/dx  is  the  reciprocal 
of  df(x)/dx,  it  follows  that  the  root  of  the  equation  x  =  f(x)  can  be 
found  either  by  the  sequence 

r-2  =  /(>'i)>     >3  =  /O  a)»     r*  =  f  0  3). 
or  by  the  sequence 


provided  only  that  /'  (x)  is  positive  and  different  from  unity  near  the  re- 
quired root. 

To  more  equations  than  one,  even  to  a  set  given  in  such  a  form  as 

x  =  f(x,  //,  z),      ij  =  g(x,  y,  z\     z  —  h(x,  y,  z), 

this  process  is  not  necessarily  adaptable,  but  the  knowledge  of  the  process 
may  well  be  an  encouragement  to  the  use  of  very  crude  approximations  in 
choosing  the  rational  coefficients  which  we  have  denoted  in  general  by 

•  ••>  Pun- 


IV. 

In  illustration  of  the  general  method  we  have  described,  we  give  some 
details  of  the  calculation  in  the  case  of  the  equations  connected  with  the 
covering  puzzle.  First,  a  rough  calculation,  of  which  no  account  need  be 
given,  shows  that  86°,  36°,  72°,  66°  is  a  better  approximation  than  36°, 
36°,  72°,  72°  to  the  solution  of  the  equations.  To  have  only  acute  angles 
to  consider,  we  substitute  36°  —  a,  36°--2/3,  72°+y,  66°—  n  for  6,  <j>,  \js,  x, 
and  the  equations  solved  are 


-2  sin  (86°—  a)    +cos(86°—  u-/8  +  y)  +  sin(180+"+^+y)  =  0, 
-2  sin  (86'  -2/3)+  cos  (30°—  a  —  /8-)7)+sin(12°+u+/8-;/)  =  0, 

-2  sin  (36    -a)  —  cos(12°—  a  —  >/)  +  2  sin(2u  +  y)  +  siii(800+u—  >?) 

+cos(18°-2/3  +  y)  +  sin(36°+2/3  +  y)  =  0, 

)  —  cos(6°  —  2/3  —  2y  +  >/)  +  cos(42°—  2u  —  2y+>/) 

=  0. 


1914.]         THK  SOLUTION  OF  NUMERICAL  FUNCTIONAL  KQUATIONH.  822 

The  equations  which  the  classical  process  requires  us  first  to  solve  are 

-3fi3y'  ='058, 

=  "102, 
'  =  '115, 
=  '02(5, 

a,  /3,  y,  i;  being  expressed  in  radians  ;  no  attempt  is  made  to  solve  these 
equations,  hut  the  set 


that  is  to  say 

19a  +   9/3+   2y  =-    Qt, 


25a+ 


is  inverted,  not  accurately,  but  to  the  approximate  form 


the  form*  of  coefficient  used  here  being  the  most  convenient.     Substituting 

-'058,  —'102,    -'115,  —'026  for  t,  u,  v,   w,   for  a  second  approxima- 

tion   to   the    values  of  a,   (3,   y,   >;  (simultaneous  zeroes   being    the    first 


*  If  a  set  of  equations  with  integral  coefficients  is  inverted  accurately,  the  coefficients  in 
the  reciprocal  set  are  fractions  with  a  common  denominator  A,  the  determinant  of  the 
original  coefficients,  and  if  the  set  of  numbers  A,  A/2,  A/3,  ...,  A/(&  —  1)  is  written  down,  k 
being  any  integer  greater  than  9,  an  approximation  of  any  desired  accuracy  to  each  coefficient 
in  the  convenient  form 


,,,  .. 

s,        10s,       100s:j       1000s4 

where  each  of  the  letters  e,,,  el}  e.<,  ...  stands  for  one  of  the  three  numbers  1,  U,  —  1,  and  each 
of  the  letters  s,,,  s,,  s;,  ...  for  a  positive  integer  less  than  k,  can  be  written  down  at  sight. 
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approximation),  we  have* 


29 


70 

4  84 

4 


29  1 

1                            17  7 
217 

26                             15  11 

a2  —      '029           /32  =  -  '004           y.2  =  —  '058  »;,  =      '058 

54  19 

15  1 

-       1  -     58 


—     21  —     73 

that  is, 

«2  =  -008  =  28',    /32  =  -021  =  1°12',    y2  =  '012  =  41',    %  =  '006  =  21'. 

To  find  a  third  approximation,  the  values  of  the  functions  on  the  left-hand 
sides  of  the  equations  for  the  values  a2,  /32,  y2>  »/2  of  the  variables  are  found, 
this  time  to  five  places  of  decimals  ;  the  calculation  requires  only  addition 
and  subtraction  of  numbers  read  from  tables,  and  the  values  found  are 
•00439,  '00663,  '00516,  --'00170.  By  substitution  of  these  values  for 
t,  u,  v,  10  in  the  formulae  already  used,  the  third  terms  in  the  approxima- 
tion are  found  to  be 

a3  =  -  -  "00094  =  -  3',         /33  =  -  '00127  =  --  4', 
73  =      -00090  =       3',  »73  =  —  '00087  =  —  3', 

and  these  are,  in  fact,  true  to  a  single  minute.  For  continuing  the 
approximation,  the  linear  equations  to  which  the  functional  equations  are, 
to  seven  places  of  decimals,  now  equivalent,  are  written  down,  the  variables 
being  expressed  in  seconds  and  denoted  by  a,  b,  c,  e,  and  the  coefficients 
being  found  from  the  difference  columns  in  the  tables.  For  example,  the 
term  cos  (36°— a— /3+y)  in  the  first  function  is  replaced  by 

cos  |360-(«,-f/9.2-y2)-(a3+^3-y3)-(*+6-c)|, 

:  This  arrangement  of  the  figures,  though  not  elegant  in  appearance,  is  the  most  con- 
venient in  practice  :  the  terms  are  calculated  in  order,  and  each  one  is  placed  above  or  below 
those  already  written  according  as  it  is  positive  or  negative.  Separate  addition  of  the  positive 
and  negative  components  is  desirable  on  account  of  the  uncertainty  of  the  sign  which  is  to 
prevail.  After  this  we  give  only  the  results  at  the  various  stages,  but  nothing  is  actually 
omitted  except  such  columns  as  these. 
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that  is,  by   cos  [35°11'  —  (<7  -\-b  —  r)\  ,    and  this  in  turn  by 
|8173125  +  (1676/60Hr/  +  /;-r);  XlO"7. 
The  equations  so  found,  multiplied  by  fiOxlO7,  are 

9137«  +   44056+1053c  -121500, 

-1459?  =  •  -  155580, 
252r  —  3085f  :       -182040, 
4610«-     55976  +  9683c  +  6666e  =  -      41460, 


equations  in  which  the  coefficients  of  the  variables,  divided  by  2909,  that 
is,  by  107  times  the  number  of  radians  in  a  minute  of  arc,  necessarily 
differ  but  little  from  the  rational  coefficients  in  the  earlier  equations  con- 
necting t,  it,  v,  w  with  a,  ft,  y,  »/.  Four  equations  such  as 


connecting  four  variables  k,  I,  m,  n  with  the  four  variables  a,  b,  c,  e  are 
therefore  approximately  equivalent  to  four  equations  of  which  the  first  is 

2909«  =  -Q-ty 


First  approximations  to  a,  b,  c,  e,  or  fourth  approximations  to  a,  ft,  y,  >/, 
are  immediately  found  to  be  -9",  -10",  1",  —10";  these  cannot  be 
trusted  to  a  second,  but  to  proceed  one  stage  further  it  is  necessary  only 
to  find  four  numbers  7i'2,  /.2,  m.2,  «2  by  substituting  these  last  terms  in  such 

equations  as 

i  +  1053^  +  121500, 


and  then  to  find  the  final  terms  in  the  approximation  from  such  equations 


It  is  found  that,  to  two  significant  figures, 

fc2  =  --  3700,     L2  =  -  6400,     m.z  =  -  6300,     wa  =  -  -  1000, 

and  that  b.2  is  slightly  less  than  ^",  while  «2>  C2>  e-z  are  about  1",  i",  J". 
The  conclusion  is  that  the  required  solution  of  the  given  equations,  correct 
to  the  nearest  second,  is 

0=  35°35'9",     «/>  =  33°  44'  19",     ^=  72°44'1",     x  =  65°  42'  10". 


The  value  of  0  +  ^0,  the  angle  whose  cosine  is  the  ratio  of  the  radii  of 
the  circles,  is  52°27'18^",  accurate  in  point  of  fact  to  at  least  ^",  and  the 
cosine  is  "6094183,  with  an  error  less  than  2  in  the  last  place  :  we  may 
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say  confidently  that  covering  is  possible  if  the  ratio  exceeds  '6094185, 
impossible  if  the  ratio  is  less  than  '6094180.  The  ratio  of  a  —  c  to  b  is 
sinl°33'43^"/sm72044'l",  that  is,  '028545. 

The  smallness  of  (a—c)/b  suggests  a  supplementary  question,  that  of 
the  least  value  of  afb  for  which  covering  is  possible  while  three  of  the 
small  circles  actually  pass  through  the  centre  of  the  large  circle.  To 
solve  this,  we  have  only  to  take  instead  of  the  equation  expressing  that 
20-J-0  is  a  maximum,  the  geometrical  condition  that  C  and  K  coincide,  a 
condition  expressed  by  the  equations 

x  =  ^  =  2<j>  =  20. 
The  angle  6  must  satisfy  the  equation 

2cosO  cosfO  =  1, 
that  is, 

1  —  cos  f  6  —  cos  %6  =  0, 

so  that  %6  satisfies  the  equation 

(g)  =  1—  cos  5^—  cos  £  =  0, 


a  first  approximation  to  the  required  root  of  this  equation  being  18°.  A 
good  approximation  to  l/(5  sin  90°  +  sin  18°)  is  (i—  T^),  and  the  root  is  the 
sum  of  a  series  fi  +  ^+.--,  whose  terms  are  found  in  succession  from 
the  formulae 


It  is  found  that  in  fact  F(g  i  +  £2  +  £3)  vanishes  to  seven  places  of  decimals, 
^1  +  ^2+^3  being  equal  to  17°28'16i",  so  that  the  ratio  of  the  radii,  being 
the  cosine  of  |0,  that  is,  of  52°24'48|",  is  '6099579. 

Perhaps  the  most  curious  feature  of  the  whole  problem  is  the  nearness 
of  this  ratio,  on  the  one  hand  to  the  smallest  ratio  permitting  covering, 
and  on  the  other  hand  to  the  smallest  ratio  which  allows  the  five  small 
circles  all  to  pass  through  the  centre  of  the  large  circle,  this  last  ratio 
being  \  sec  36°,  that  is,  '6180340.  The  difference  between  the  smallest 
ratio  allowing  covering  and  the  last  ratio  found  is  just  large  enough  to 
take  effect  in  practice,  but  the  lack  of  precision  in  the  painted  circle  and 
the  thickness  of  the  discs  prevent  the  accuracy  which  would  be  necessary 
if  a  distinction  was  to  be  made  between  the  first  two  arrangements  dis- 
cussed, and  if  my  readers  can  perceive  the  centre  of  the  large  circle  they 
may  proceed  to  pocket  as  many  of  the  showman's  rewards  as  they  feel 
themselves  to  have  earned  by  reading  these  pages. 

[My  brother,  Mr.  B.  M.  Neville,  to  whom  I  am  indebted  for  the  drawing 


1914.]          THE  SOLUTION  OF  NUMERICAL  FUNCTIONAL  EQUATIONS. 


82U 


to  scale  which  accompanies  this  paper,  finds  that  unless  the  large  circle 
has  a  diameter  of  about  a  metre  the  possibility  of  the  best  construction, 
when  the  construction  through  the  centre  is  inadequate,  cannot  be  made 
evident.  To  him  is  due  the  discovery  of  a  most  convenient  and  accessible 
form  of  covering  disc— the  pieces  of  parchment  sold  as  jam  covers  :  pro- 
vided with  five  of  these,  one  has  only  to  draw  on  a  sheet  of  paper  a  large 
circle  of  appropriate  radius  ;  the  transparency  of  the  small  circles  is  a 
great  advantage.] 
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